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1.  INTRODUCTION 


The  main  goal  of  this  research  effort  is  to  establish  a  master  dose  response  curve  correlating  blast 
overpressure  with  mortality  rate  using  a  field  validated  blast  injury  animal  model.  This  curve  will 
map  the  probability  of  24  hour  survival  over  a  wide  range  of  field-relevant  blast  over  pressures.  We 
used  peak  overpressure  (0-450  kPa)  and  impulse  (0-800  Pa-s)  as  two  predictors  brain  injury,  and 
based  on  these  values  established  the  thresholds  for  blast  TBI  in  rat  model.  This  model  was  further 
correlated  with  systemic  pathophysiological  changes:  bradycardia  and  extent  of  lung  injury,  to 
determine  if  these  changes  can  be  used  as  convenient  markers  of  injury  which  can  be  used  in  the 
field  and  inform  affected  personnel  immediately  upon  the  exposure  to  shock  wave.  We  used 
experimental  measures  to  probe  the  response  of  the  body  to  blast  loading  via  implantation  of 
intracranial  pressure  (ICP)  and  carotid  artery  (CA)  sensors  in  anesthetized  live  animal  model.  The 
relationship  between  the  blast  over  pressures  (external  loading  to  the  head)  and  the  resulting 
waveforms  of  intracranial  and  arterial  pressure  (internal  loading  to  the  brain  and  torso)  indicates 
wave  is  traveling  in  the  body  gradually  losing  its  potential.  We  have  initiate  the  development  of 
Rat  Testing  Device  (RTD),  a  prototype  surrogate  model  equipped  with  variety  of  pressure  sensors 
and  accelerator  to  precisely  determine  loading  conditions  and  facilitate  comparison  in  different 
blast  loading  scenarios.  Developed  numerical  model  of  the  shock  wave  propagation  in  the  shock 
tube  indicates  excellent  correlation  with  experimental  data.  Numerical  model  of  the  rat  helped 
identify  loading  pathways  through  which  blast  overpressure  ‘leaks’  into  the  brain  leading  to 
differential  loading  of  the  tissue  in  various  discrete  regions. 


2.  KEYWORDS 


3.  ACCOMPLISHMENTS 


Major  Goals  of  the  Project  (Statement  of  Work  with  Timeline): 


1  STATEMENT  OF  WORK 

Year  1 

Year  2 

Year  3 

Year  4  [ 

Project  Steps 

|Q1|Q2|Q3|Q4 

|Q'|Q2|Q3|Q4 

I  W  1  Q2 1  Q.q  Q-t 

Wlw|Q3|Q4| 

SA  1:  Develop  Master  Dose  Response  Curves  for  10  week  old  SD  rats 

Task  1:  Evaluate  Mortality  Rates  and  Biomechanical  Loading  in  Wide  Range  of  Blast  Intensities 
Task  2:  Determine  the  Biomechanical  Loading  of  the  Rat  Brain  During  Simulated  Blast 
Task  3:  Numerical  Simulation  of  Brain  Injury 
Milestone  1:  Master  Dose  Response  Curve  Developed 
SA  2:  Assess  Pathologies  in  Mild-to-Moderate  bTBI  Range  24  hours  After  the  Exposure 
Task  4:  Evaluate  the  Lung  Injuries  Caused  by  Blast  Exposure 

Task  5:  Assess  Extent  of  Oxidative/Nitrosative  Stress,  BBB  Damage  and  Neuroinflammation 
Task  6:  Examine  Plasma  Membrane  Damage  Using  Fluorescent  Tracers 
Task  7:  Evaluate  Alterations  in  Brain  Proteome  After  Primary  Blast  Exposure 
Milestone  2:  Brain  Injury/Blast  Load  Correlated  with  Bioniarkers 
SA  3:  Examine  the  Effect  of  Blast  Impulse  on  Master  Dose-Response  Curve 

Task  8:  Establish  Master  Impulse  Dose-Response  Curve  at  Three  Blast  Overpressures 
Task  9:  Examine  Changes  in  Protein  Expression  Due  to  Changes  in  Blast  Impulse 
Task  10:  Determine  Alterations  of  Loading  in  the  Rat  Brain  Caused  by  Changes  in  Impulse 
Milestone  3:  The  Effects  Associated  With  Blast  Impulse  Identified 
SA  4:  Establish  Human  Injury  Criterion  for  Blast  TBI  -  the  Cross-Species  Correlation  Function 
Task  1 1 :  Identify  Cross-Species  Milt-To-Moderate  bTBI  via  Numerical  Simulations 
Task  12:  Comprehensive  Characterization  of  Blast  Loading  in  PMHS  Head  Model 
Task  13:  Devise  Transfer  Function  to  Correlate  Rat  and  Human  Blast  Overpressure  Conditions 
Milestone  4:  The  Human  Injury  Criterion  Established 
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Completion:  Task  1  and  2  completed;  Task  3  is  in  progress  (50%  complete),  Tasks  4  &  5  (25% 
complete) 


Task  1:  Evaluate  Mortality  Rates  and  Biomechanical  Loading  in  Wide  Range  of  Blast  Intensities 


Optimization  of  9-inch  shock  tube  operations 

The  comprehensive  calibration  of  the  9-inches  square  cross  section  shock  tube  was  completed, 
including  data  analysis  and  manuscript  preparation.  The  new  9-inch  square  cross  section  shock 
tube  was  designed  and  fabricated  based  on  extensive  testing  experience  from  the  past,  but  there 
is  still  room  for  improvement.  For  example,  during  recently  performed  experiments  we  have 
optimized  the  pressure  sensor  mounting,  which  otherwise  might  sometimes  introduce  artifacts 
into  the  recorded  pressure  history.  Armed  with  this  new  experience  we  are  constantly  improving 
our  pressure  measurement  procedures  to  get  as  accurate  measurements  as  possible. 

Additionally,  the  modifications  introduced  during  manufacturing  process  as  well  as 
environmental  factors  which  might  influence  the  profiles  of  generated  shock  waves  compared  to 
original  design  and  it  is  more  than  likely  due  to  manufacturing  process  some  parts  critical  for 
shock  tube  operation  are  different,  which  might  influence  its  performance.  Thus  the 
comprehensive  calibration  gave  us  an  insight  into  the  performance  of  the  shock  tube  and  helped 
establish  a  basic  set  of  working  parameters.  In  the  same  time  there  were  a  number  of  parameters 
and  new  sensors  which  were  tested  and  optimized  in  order  to  make  these  measurements  as 
reliable  as  possible.  For  this  purpose  we  performed  additional  test  shots  to  answer  questions 
regarding  the  source  of  observed  overpressure  noise  and  aberrations  adversely  affecting  recorded 
waveforms. 

Dose-response  logistic  regression  model 

All  the  test  animals  were  in  a  prone  position  with  a  head  restrained  and  thus  no  artifacts 
associated  with  uncontrolled  head  acceleration  were  present  in  this  bTBI  model,  as  demonstrated 
elsewhere.1  The  net  movement  of  the  specimen  did  not  exceed  3  mm  as  determined  by  analysis 
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of  captured  high  speed  video.  We  exposed  13  rats  per  discrete  incident  peak  overpressure  to 
develop  a  dose-response  linear  regression  model  for  predicted  mortality  rate  (PMR)  as  a  function 
of  BOP  and  corresponding  impulse  values  (fig.  1).  We  did  not  observe  any  mortality  of  animals 
exposed  to  blast  intensity  lower  than  170  kPa  (fig.  1,  left)  or  impulse  300  Pa-s  (fig.  1,  right). 
However  mortality  rate  gradually  increased  from  170  kPa  to  300  kPa,  and  mortality  rate  was 
100%  above  300  kPa  (or  500  Pa-s  impulse).  These  accounted  for  an  immediate  loss  (not  delayed 
deaths,  i.e.  animals  were  removed  from  the  shock  tube  in  less  than  a  minute  after  the  shock  wave 
exposure)  of  45  rats  as  a  consequence  of  shock  wave  exposure,  consistent  with  our  previous 
work.  ’  The  insets  in  fig.  1  present  receiver  operating  characteristic  (ROC)  curves  for  respective 
logistic  regression  fits  with  the  following  areas  under  ROC  curves:  0.878  (left)  and  0.882  (right), 
for  BOP  and  impulse,  respectively.  Analysis  from  these  models  indicated  good  predictive  power 
of  mortality  rate  using  peak  overpressure  and  impulse  as  metrics  to  gauge  injury  risk  under 
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primary  blast.  The  McFadden  pseudo-R  values  of  0.389  and  0.412  were  obtained  for  fits  using 
peak  overpressure  and  impulse  as  independent  variables,  respectively.  These  values  indicate  a 
highly  satisfactory  quality  of  fit. 

The  logistic  regression  models  were  developed  and  evaluated  using  Systat  13.0  software  (Systat 
Software,  Inc.,  San  Jose,  CA).  Dose-response  models  for  heart  rate  and  pulmonary  injury  were 
fitted  with  Origin  9.0  software  (OriginLab  Corp.,  Northampton,  MA)  using  dose-response 
function: 

y  =  A±  +  (A2  -  AJ/il  +  10(logx°-*)P)  (1) 

Where  Aj,  A2  are  asymptotes,  log  x0  is  an  inflection  point,  and  p  is  a  slope  value. 


Figure  1  The  logistic  regression  dose-response  model  for  rats  exposed  to  single  blast  with  intensity  in  the  range 
of  60-450  kPa  peak  overpressure  (left)  and  corresponding  impulse  (90-780  Pa-s,  right).  The  insets  present 
receiver  operating  characteristic  (ROC)  curves  for  respective  logistic  regression  fits  with  areas  under  ROC 
curves:  0.88  and  0.88.  McFadden  R2  =  0.39  and  0.41,  for  fits  using  peak  overpressure  and  impulse  as 
independent  variables,  respectively.  Dashed  lines  indicate  confidence  intervals  of  the  model. 


Primary  blast  impacts  mild  pulmonary  injury 

Animals  tested  in  prone  position  have  their  abdomen  partially  protected  from  the  blast  wave  by 
the  aluminum  holder  used  in  our  experiments.1  The  observed  levels  of  pulmonary  injury 
expressed  using  Yelveton’s  scoring  systenr  revealed  a  low  level  of  injury  (fig.  2  A,  B). 
However,  we  observed  an  increasing  trend  of  injury  score  with  increasing  peak  overpressure  and 
impulse.  We  observed  only  a  few  cases  where  pathological  score  exceeded  21  for  the  blast 
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strength  higher  than  300  kPa  BOP  with  high  standard  deviations  (fig.  2  C).  A  score  of  21  is 
considered  as  a  cut-off  threshold  for  mild  pulmonary  injury.3  The  pathological  score  at  50% 
PMR  (at  260  kPa)  was  found  to  be  less  than  10,  while  the  score  was  less  than  4  in  the  60  -  190 
kPa  range.  Moreover,  there  are  six  animals  which  died  as  a  result  of  blast  exposure  and  had  no 
lung  injury  (score  of  zero,  fig.  2C,  D).  These  results  suggested  minimal  pulmonary  injuries  and 
thus,  we  conclude  lung  injury  is  not  a  viable  indicator  of  PMR. 


Blast  induced  bradycardia 


We  evaluated  the  functional  changes  in  the  heart  rate,  blood  oxygen  saturation  (spCL)  and 
perfusion  index  in  the  60-250  kPa  peak  overpressure  range  over  period  of  30  minutes  before  and 
after  blast  exposure.  We  found  that  the  onset  of  bradycardia  occurred  immediately  after  the  blast 
exposure  even  at  60  kPa  (fig.  4,  p=0.01,  power:  0.85  vs  control).  The  values  of  the  differential 
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Figure  2  Lung  injury  scores  for  rats  exposed  to  a  single  blast.  The  dose-response  model  was  used  to  fit  the  IS  as 
a  function  of  peak  overpressure  (A,  C)  and  impulse  (B,  D).  The  scattergrams  (C,  D)  illustrate  individual  scores 
and  their  distribution  among  the  cohort  of  75  rats  evaluated  in  this  test.  The  value  of  21  is  the  upper  limit  of  the 
slight  lung  injury  level  as  defined  by  Yelverton.  The  vertical  lines  (peak  overpressure  of  260  kPa  (C),  or  impulse 
of  450  Pa-s  (D))  correspond  to  50%  predicted  mortality  rate  according  to  the  dose-response  model  in  fig.  1. 
There  are  six  animals  which  died  after  the  blast  and  had  no  apparent  lung  injuries  (score  of  zero). 


average  heart  rates  (AHR)  decreased  gradually  with  increase  in  blast  intensity:  -29+10  (60  kPa),  - 
26+20  (100  kPa),  -43+26  (130  kPa),  -62+21  (190  kPa),  -62+43  (230  kPa)  and  -62+24  (250  kPa) 
bpm.  These  data  were  modeled  using  a  simple  dose-response  function  to  quantify  the 
characteristics  of  AHR  as  a  function  of  two  blast  parameters,  the  peak  overpressure  and  the 
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impulse.  This  mathematical  modeling  generated  asymptotic  values  of  Ai  and  A2:  -13.9  and  -59.5 
(for  peak  overpressure,  fig.  2A),  and  -12.7  and  -67.1  (for  impulse,  fig.  2B),  respectively.  The 
calculated  inflection  points  (log(vo))  are  109.3  kPa  and  186.9  Pa-  s  for  both  blast  parameters  (fig. 
3).  The  control  group  was  not  correlated  with  any  of  the  blast  exposure  groups  (p  <  0.05). 


Peak  overpressure  (kPa)  Impulse  (Pa  s) 

Figure  3  Blast  induced  bradycardia  as  a  function  of  peak  overpressure  (A)  and  impulse  (B).  In  both  cases,  the 
dose-response  function  (1)  was  used  to  model  the  pathological  response  and  resulting  parameters  are  listed  in 
respective  insets  in  both  figures.  All  blast  exposed  groups  have  statistically  significant  heart  rate  decrease  versus 
control  (p  <  0.05). 


Induction  of  oxidative/nitrosative  stress  markers 

Using  our  logistic  regression  risk  injury  model,  we  defined  the  upper  level  of  mTBI  at  145  kPa 


Figure  4  Mild  TBI  range  of  blast-wave  exposure  induces  NADPH  oxidase  1  (NOX1),  inducible  nitric  oxide 
synthase  (iNOS),  4-hydroxynonenal  (4HNE)  and  3-nitrotyrosine  (3NT)  expression  in  rat  brain  microvessels.  (A) 
A  representative  of  immunofluorescent  staining  of  specific  marker  as  indicated  on  the  Y-axis  in  intact 
microvessels  of  brain  cross  sections  from  rats  subjected  to  a  single  exposure  to  60,  100,  or  130  kPa  peak 
overpressure,  and  control.  (B)  Corresponding  Western  blot  specific  marker  (as  indicated  on  the  Y-axis)  and 
housekeeping  protein,  P-actin.  (C)  The  quantification  of  the  results  of  the  immunoreactive  fluorescence 
intensities.  Values  are  mean  +  SEM  (n  =  4)  with  p-value  <  0.01  compared  with  control. 
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BOP  as  5%  PMR.  We  examined  the  oxidative/nitrosative  injury  and  the  BBB  integrity  in  the 
range  of  60  kPa  to  130  kPa  peak  overpressures  and  determined  lower  threshold  of  mTBI  at  80 
kPa.  These  markers  were  evaluated  24  hours  after  the  injury,  at  a  time  point  which  proved  in  our 
previous  study  to  yield  their  maximum  levels  in  analogous  bTBI  model.4  We  first  evaluated  the 
induction  of  free  radical  generating  enzymes  NADPH  oxidase  (NOX1)  and  inducible  nitric  oxide 
synthase  (iNOS)  in  the  brain  capillary  cross  section  tissues.  We  found  that  blast-wave  exposure 
significantly  up-regulated  the  induction  of  NOX1  with  100  kPa  (p=  0.02)  and  130  kPa  (p=  0.001) 
BOP  (fig.  4,  left  upper  panel).  Similarly,  iNOS  expression  was  increased  at  130  kPa  BOP  (p= 
0.0005)  and  at  100  kPa  BOP  (p=0.04)  (fig.  4,  left  lower  panel).  Induction  of  NOX1  and  iNOS 
produces  superoxide  and  nitric  oxide  respectively,  which  will  also  react  together  to  form 
peroxynitrite,  a  more  reactive  free  radical.  The  oxidative/nitrosative  damage  is  a  post-oxidant 
production  event.  Proteins  adducted  with  4-hydroxynonenal  (4HNE,  oxidative  stress  marker)  or 
3-nitrotyrosine  (3NT,  nitrosative  stress  marker)  are  used  for  assessing  the  extent  of  injury  in  the 
tissue.  In  parallel  with  the  induction  of  NOX1  and  iNOS  enzymes,  we  found  that  the  level  of 
oxidative  damage  signature  4HNE  (fig.  4,  right  upper  panel)  is  increased  in  mTBI  range 
exposure  with  100  kPa  (p=  0.05)  and  130  kPa  (p=  0.0004).  Similarly,  3NT  (fig.  4,  right  lower 
panel)  was  found  to  be  increased  with  the  mTBI  range  exposure  (100  kPa  BOP,  p=  0.04  and  130 
kPa  BOP,  p=  0.002). 


Disruption  of  the  BBB  integrity 


The  capillary  oxidative/nitrosative  damage  might  lead  to  BBB  disruption,  and  thus  we  evaluated 
the  alterations  of  tight  junction  (TJ)  proteins  claudin-5,  occludin  and  zonula  occluden  1  (ZO-1). 
TJ  proteins  are  the  primary  functional  barrier  biomolecules  of  the  BBB.  A  reduction  in  TJ 
protein  levels  or  disruption  of  the  architectural  structure  of  any  TJ  protein  is  expected  to  impair 
BBB  integrity,  thereby  enhances  the  chance  of  immune  cell  infiltration  into  the  brain  for 
initiation  of  neuroinflammation. 


Control  60  kPa  100  kPa  130  kPa 


Control  60  kPa  100  kPi  130  k Pa 


Figure  5  Oxidative/nitrosative  damage  of  capillaries  impaired  BBB  integrity  by  disrupting  the  tight  junction 
proteins.  (A).  A  representative  of  immunofluorescent  staining  of  TJ  protein  Claudin-5  (a),  Occludin  (b),  and 
zonula  occluden  1  (ZO-1,  c)  in  capillaries  of  whole  brain  cross-sections  in  control  and  blast  exposed  animals. 
(B).  Corresponding  Western  blot  and  bar  graphs  showing  the  quantification  of  respective  immunoreactive 
fluorescence  intensity  of  Claudin-5  (d,  e),  Occludin  (f,  g),  and  zonula  occluden  1  (h,  i).  Values  are  mean  ±  SEM, 
(n  =  4),  and  asterisk  indicates  each  blast  intensity  compared  with  control  where  level  of  statistical  significance 
was  achieved  (p  <  0.05). 
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Interestingly,  our  results  showed  that  mTBI  range  of  blast-wave  exposure  decreased  the  levels  of 
claudin-5  (130  kPa,  p  =  0.001)  (fig.  5  a,  d,  e),  occludin  (130  kPa,  p=  0.002)  (fig.  5  b,  f,  g)  and 
ZO-1  (100  kPa,  p  =  0.03;  130  kPa,  p  =  0.0007)  (fig.  5  c,  h,  i)  proteins  in  brain  tissue  sections. 
These  data  suggest  possible  leakiness  of  the  BBB  and  it  might  result  in  neuronal  inflammation 
around  the  perivascular  region  of  the  brain. 

Assessment  of  BBB  leakage 

To  assess  this  BBB  leakiness  and 
neuronal  injury,  we  examined  the 
leaking  out  of  neuronal- specific 
enolase  (NSE)  into  the  blood 
samples  of  tissues  exposed  at  BOP 
of  60,  100  and  130  kPa  versus 
control  animals.  In  agreement  with 
a  decrease  in  BBB  integrity,  we 
observed  elevation  of  NSE  levels 
in  blood  samples  exposed  to  blast 
compared  with  controls  (fig.  6). 

The  glycolytic  enzyme  enolase  is  a 
dimeric  isoenzymes  and  also 
known  as  neuronspecific  enolase 
(NSE-aa,  ag  and  gg),  as  these 
isoenzymes  were  initially  detected 
in  neurons  and  neuronendocrine 
cells.  However,  in  other 
pathological  conditions  such  as 
small  cell  lung  cancer  and 
neuroblastoma,  NSE  exhibits  a 
signatory  value  in  disease 
detection  and  progression.  In  the 
absence  of  such  pathological  status  or  stimuli,  as  in  our  experimental  setup,  NSE  detection  in 
plasma  samples  establishes  BBB  leakage  in  primary  bTBI. 

Aperio  VERSA  200  slide  scanner 

The  installation  of  the  new  casework  in  the  biochemistry  suite  was  completed  (under  the 
financial  support  from  internal  funds  of  NJIT).  We  thus  proceeded  with  the  installation  of  the 
Aperio  VERSA  200  slide  scanner  (fig.  7).  The  installation  of  the  hardware  and  related  training 
sessions  were  performed  in  March  2016.  Subsequently  we  have  acquired  the  server  to  run  the 
database  and  analysis  eSlide  Manager.  The  system  was  ensured  to  be  fully  operational  and  tested 
before  advancing  to  the  next  step,  i.e.  training  program  utilizing  two  acquired  analysis 
algorithms:  1)  Area  Quantification  FL,  and  2)  Genie.  These  activities  were  concurrent  with 
scanning  of  slides  prepared  by  members  of  the  group  to  familiarize  themselves  with  the  system 
in  order  to  perform  planned  research  tasks  efficiently.  To  date  more  than  200  slides  with  coronal 
sections  were  digitized  using  this  system  and  it  will  be  the  core  of  our  activities  in  the  future. 


Figure  6  The  neuronal  damage  marker  leakage  into  serum  is 
facilitated  by  BBB  disruption.  ELISA  results  show  the  levels  of 
neuronal  specific  enolase  (NSE)  in  rat  blood  serum  collected  from  a 
single  blast  exposure  at  60,  100  and  130  kPa.  Blood  samples  were 
collected  at  the  time  of  sacrifice  after  24  hour  post  blast  exposure. 
Values  are  mean  +  SEM,  n=4.  Statistical  significance  (p  <  0.01) 
compared  with  control  was  observed  at  130  kPa  as  the  lowest  blast 
intensity. 
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Figure  7  The  overview  of  the  new  Aperio  VERSA  200  slide  scanner  system:  A)  Dell  PowerEdge  T630  server,  B) 
the  workstation  and  NAS  used  to  collect  the  data,  and  C)  the  slide  scanner  performing  data  acquisition.  The 
stained  and  digitized  coronal  section  of  the  rat  brain  is  visible  on  the  monitor. 


Task  2:  Determine  the  Biomechanical  Loading  of  the  Rat  Brain  During  Simulated  Blast 

Correlation  of  intracranial  pressure  during  blast  exposure  and  brain  injury  mechanisms  is 
ongoing  research  topic  where  many  important  questions  are  awaiting  to  be  appropriately 
addressed.  Particularly,  due  to  limited  number  of  studies  very  little  experimental  evidence  is 
available  in  the  literature  correlating  specific  brain  pathologies  with  ICP  fluctuations.  To  date 
most  of  the  experimental  work  was  performed  using  numerical  models  and  direct  ICP  pressure 
measurements.  These  studies  give  only  limited  insight  into  BINT  associated  with  ICP 
fluctuations.  The  contribution  to  brain  injury  associated  exclusively  with  skull  flexures  is  a 
difficult  research  task  aggravated  by  simultaneous  occurrence  of  other  injury  factors  contributing 
to  increase  of  ICP  during  blast  exposure. 

Experiments  performed  in  our  laboratory  at  University  of  Nebraska-Lincoln  preformed  on  10 
weeks  old  Sprague  Dawley  rats  exposed  to  five  discrete  blast  intensities  (130,  190,  220,  250  and 
290  kPa  peak  overpressure)  resulted  in  peak  incident-to-ICP  ratios  in  the  range  of  1.5  to  1.3  for 
lower  and  higher  peak  overpressures  (see:  Skotak  et.  al.,  J.  Neurotrauma,  2013,  30,  1147-1160). 
The  FFT  analysis  of  ICP  profiles  indicated  the  oscillations  have  harmonic  characteristics,  with 
low  (5-10  kHz)  and  high  frequency  component  (10-20  kHz).  The  high  frequency  oscillations 
were  linked  to  observed  mortality  (animals  exposed  to  190,  250  and  290  kPa  BOP)  and  we  have 
also  observed  the  maximum  ICP  frequency  increases  with  increasing  blast  intensity:  13  kHz  (190 
kPa),  16  kHz  (250  kPa)  and  19  kHz  (290  kPa). 


Time  (ms)  Nominal  pressure  Nominal  pressure 

Figure  8  Representative  pressure  profiles  measured  by  incident  (loading),  intracranial  and  carotid  artery  sensors 
(response)  implanted  in  a  rat  exposed  to  130  kPa  of  nominal  blast  intensity.  Biomechanical  loading  was 
performed  at  5  blast  overpressures  in  the  100  to  250  kPa  range  and  peak  overpressure  (middle  panel)  and 
impulse  (right  panel)  were  quantified. 
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However,  our  measurements  performed  on  a  cohort  of  6  animals  indicate  ICP-to-incident 
pressure  ratio  is  less  than  1,  which  is  a  strong  contradiction  of  what  has  been  reported  in  the 
literature  to  date  (fig.  8).  Most  likely  the  methodology  differences  are  responsible  for  observed 
discrepancies.  We  have  initiated  use  of  miniature  Millar  SPR-671  sensors  which  are  only  0.5  mm 
OD  and  typically  are  suitable  for  measurements  of  arterial  blood  pressure.  These  sensors  have 
robust  construction  and  easily  withstand  blast  exposure,  which  was  confirmed  during  our 
experiments.  Currently  we  are  performing  implantation  of  sensors  using  live  animals  and  this 
might  be  a  contributing  factor  responsible  for  observed  differences.  The  sensor  size  and 
mounting  strategy  might  also  play  important  role.  In  the  past  the  sensor  was  mounted  in  the 
forehead,  which  is  not  ideal  location  considering  biofidelity  of  the  measurement,  i.e.  when  the 
skull  structure  is  affected  by  craniotomy  and  application  of  the  cement  to  keep  sensor  cannula  in 
place,  one  would  expect  the  pressure  transmission  of  the  head  to  be  quite  different  than  when 
skull  and  skin  layer  remain  intact.  It  is  also  hard  to  avoid  relative  motion  of  the  sensor  with 
respect  to  skull  using  this  strategy,  although  it  might  be  a  general  drawback  of  any  sensor 
implantation  methodology  due  to  specificity  of  blast  experimentation.  There  is  always  a  minute 
head  displacement  independently  how  well  the  restrain  system  is  designed  and  performed. 

Task  3:  Numerical  Simulation  of  Brain  Injury 
Rat  head  and  shock  tube  models 

High  resolution  MRI  and  CT  datasets  of  a  male  Sprague  Dawley  rat  were  used  in  combination  to 
generate  a  three  dimensional  rat  head  model.  Separate  T2-weighted  MRI  scans  for  muscle  and 
skin,  and  brain,  along  with  a  CT  scan  for  skull  and  bones  were  used  to  achieve  proper  contrast 
and  segmentation  of  various  tissues  (i.e.  muscle,  skin,  brain,  skull  and  bones).  Brain  MRI 
consisted  of  an  isotropic  resolution  of  256  x  256  x  256  pixels,  for  a  field  of  view  of  30  mm  in  all 
three  directions.  MRI  for  muscle  and  skin  has  an  anisotropic  resolution  of  512  x  512  x  256,  for  a 
field  of  view  of  30,  30  and  50  mm  respectively.  The  three  datasets  where  overlapped,  registered, 
segmented  and  triangulated  using  Avzio  6.2.  The  triangulated  surface  mesh  was  imported  to 
Hypermesh,  and  a  volume  mesh  consisting  of  10-node  tetrahedron  element  was  generated  from 
the  imported  mesh.  The  generated  interface  is  shared  between  the  brain,  skin  and  skull  regions, 
which  are  treated  as  Lagrangian  elements.  The  model  is  then  imported  to  Abaqus  6.13,  and 
inserted  into  the  shock  tube  model.  The  generation  and  propagation  of  blast  waves  are  modeled 
in  the  shock  tube  environment.  The  air  inside  the  tube,  in  which  the  blast  wave  propagates,  is 
modeled  with  Eulerian  elements.  The  size  of  the  tube  corresponds  to  the  physical  dimensions  of 
the  shock  tube  used  in  experiments,  with  the  cross  section  229  x  229  mm.  The  tube  was 
subjected  to  biased  meshing  (mesh  size  varying  from  100mm  to  5mm),  with  the  highest  density 
of  fine  meshes  in  the  region  encompassing  the  rat  head  (the  test  section)  and  coarse  mesh 
elsewhere,  to  reduce  the  total  number  of  elements  in  the  model,  without  sacrificing  accuracy. 

Loading,  Interface  and  Boundary  Conditions 

The  rat  head  model  is  subjected  to  the  blast  wave  from  the  frontal  direction.  A  complete  model 
of  the  shock  tube,  with  the  driver,  transition  and  extension  section,  reproducing  the  burst, 
expansion  and  development  of  the  blast  wave  would  be  time  consuming  and  computationally 
taxing,  thus  less  efficient  for  this  research  task.  Hence,  a  partial  model  with  experimentally 
measured  pressure-time  profile  is  used  as  the  pressure  boundary  condition  as  well  as  the  input, 
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with  the  numerical  model  including  only  downstream  flow  field  including  test  specimen  section 
of  the  shock  tube. 


Rut  Head  Assembly 


Surrounding’  medium 
modelled  with  Eulerian 
elements 


Lagrangian  body 
immersed  in  Eulerian 
domain 


Figure  9  Sequence  of  finite  element  modelling  is  presented:  1)  MRI/CT  scans  are  overlapped,  registered, 
segmented  and  triangulated  using  Avizo  6.0,  2)  triangulated  surface  mesh  is  imported  to  Hypermesh  to  generate 
3D  mesh  with  10-node  tetrahedron  Lagrangian  elements,  and  3)  model  is  then  imported  to  Abaqus  and 
assembled  with  the  Eulerian  shock  tube.  Numerical  boundary  conditions  are  set  for  the  shock  tube  along  with 
displacement  constrained  within  the  axis,  and  rat  model  with  displacement  constrained  in  all  three  linear 
directions  (x,  y  and  z). 


Except  for  the  face  on  which  the  pressure  is  acting,  velocity  perpendicular  to  all  other  remaining 
faces  of  the  shock  tube  is  kept  at  zero,  to  avoid  leaking  of  air  from  the  shock  tube.  This 
configuration  also  maintains  a  planar  shock  front  travelling  in  the  longitudinal  direction  with  no 
lateral  flow.  The  interface  between  all  components  (skin,  skull  and  brain)  of  the  rat  head  was 
modelled  as  tied  (i.e.  no  tangential  and  sliding)  contact.  The  nodes  on  the  bottom  and  rear  faces 
of  the  rat  head  is  constrained  in  all  degrees  of  freedom. 
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An  enhanced  immersed  boundary  method  is  used  to  provide  the  coupling  between  the  Eulerian 
and  the  Lagrangian  domains.  Here,  the  Lagrangian  region  resides  fully  or  partially  within  the 
Eulerian  region  and  provides  no-flow  boundary  conditions  to  the  fluid  in  the  direction  normal  to 
the  local  surface.  Further,  the  Eulerian  region  provides  the  pressure  boundary  conditions  to  the 
Lagrangian  region.  Thus,  a  combination  of  fixed  Eulerian  mesh  and  solid-fluid  interface 
modeling  through  the  enhanced  immersed  boundary  method  allows  for  the  concurrent 
simulations  of  the  formation  and  propagation  of  a  primary  blast  wave  in  a  fluid  medium  and 
accounts  for  the  effects  of  both  fluid-structure  interaction  and  structural  deformations  once  the 
blast  wave  encounters  a  solid.  The  interactions  (contact  conditions)  between  Eulerian  (containing 
air  and  a  propagating  blast  wave)  and  Lagrangian  regions  are  defined  using  ‘general  contact’ 
feature  (card)  in  Abaqus.  In  general,  contact  constraints  are  enforced  through  the  penalty  method 
with  finite  sliding  contact  formulation.  Various  contact  property  models  are  available  in  general 
contact.  In  the  present  work,  frictionless  tangential  sliding  with  hard  contact  is  used  as  the 
contact  property  model. 


Solution  Scheme 


The  finite  element  model  is  solved  using  the  nonlinear  transient  dynamic  procedure  with  the 
Euler-Lagrangian  coupling  method  (Abaqus  6.13).  In  this  procedure,  the  governing  partial 
differential  equations  for  the  conservation  of  momentum,  mass  and  energy  along  with  the 
material  constitutive  equations  and  the  equations  defining  the  initial  and  boundary  conditions  are 
solved  simultaneously. 

Conservation  of  mass  (continuity  equation): 

pd^+di+v^p  =  Q  ^ 

Conservation  of  momentum  (equation  of  motion): 

+  pbi  —  pat  (1.6) 


Conservation  of  energy  (energy  equation): 


de  _  dVi  dq< 

p-+We  =  oil  —  -  —  +  pqs 


(1.7) 


where,  p  is  density,  x,  v  and  a  are  displacement,  velocity  and  acceleration  of  a  particle 
respectively,  o  is  Cauchy  stress,  b  is  body  force,  e  is  internal  energy  per  unit  mass,  q  is  heat  flow 
per  unit  area  and  qs  is  rate  of  heat  input  per  unit  mass  by  external  sources. 

In  Eulerian-Lagrangian  method,  the  whole  model  is  solved  (i.e.  both  Eulerian  and  Lagrangian 
domains)  with  the  same  Lagrangian  equations.  The  notion  of  a  material  (solid  or  fluid)  is 
introduced  when  specific  constitutive  assumptions  are  made.  The  choice  of  a  constitutive  law  for 
a  solid  or  a  fluid  reduces  the  equation  of  motion  appropriately  (e.g.,  compressible  Navier-Stokes 
equation,  Euler  equations  etc.).  For  the  Eulerian  domain  in  the  model  the  results  are  simply 
mapped  back  to  the  original  mesh  with  extensions  to  allow  multiple  materials  and  to  support  the 
Eulerian  transport  phase  for  Eulerian  elements.  Eulerian  framework  allows  for  the  modeling  of 
highly  dynamic  events  (e.g.  shock)  which  would  otherwise  induce  heavy  mesh  distortion.  In 
Abaqus,  the  Eulerian  time  increment  algorithm  is  based  on  an  operator  split  of  the  governing 
equations,  resulting  in  a  traditional  Lagrangian  phase  followed  by  an  Eulerian,  or  transport 
phase.  This  formulation  is  known  as  “Lagrange-plus-remap.”  During  the  Lagrangian  phase  of  the 
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time  increment  nodes  are  assumed  to  be  temporarily  fixed  within  the  material,  and  elements 
deform  with  the  material. 

Model  validation 
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Figure  10  A.  Input  value  (incident  pressure),  B.  Unfiltered  result  at  the  hippocampus,  C.  Filtered  using 
Butterworth  filter  at  2.5  kHz,  and  D.  Comparison  of  filtered  and  unfiltered  values. 


Results  of  simulations  will  require  further  refinement,  since  excessive  oscillations  were  observed 
within  the  pressure  profiles  obtained  from  the  hippocampus  region  (fig.  10).  The  existence  of  the 
oscillations  may  be  due  to  variation  between  the  geometries  of  the  Lagrangian  vs  the  Eulerian 
domain,  leading  to  inaccurate  transfer  of  the  blast  wave  from  the  air  into  the  rat  head.  In  order  to 
eliminate  the  oscillations,  the  Butterworth  filter  was  applied,  which  led  to  the  peak  pressure 
falling  by  half.  Variation  of  mesh  size  (5  mm  to  1  mm),  made  the  oscillations  higher.  Also, 
literature  was  used  to  obtain  alternate  material  properties  for  the  skin,  skull  and  brain,  which  was 
used  for  simulations.  These  modifications  had  no  effect  on  the  results  of  repeated  simulations. 


Wave  transmission  pathways 

Numerical  model  allowed  us  to  identify  loading  pathways  in  the  rat  brain  during  the  blast 
exposure.  As  the  blast  wave  impinges  the  rat,  the  wave  first  interacts  with  the  snout  and 
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undergoes  diffraction,  where  it  bends  and  converges  towards  the  eye  socket  (pathway  1)  and  on 
top  of  the  skull  (pathway  2).  The  surface  pressure  loading  along  pathways  1  and  2  are 
transmitted  to  the  rat  brain,  which  start  converging  towards  each  other  in  the  region  of  bregma, 
lambda  and  midline  sutures.  The  loading  through  the  snout  (pathway  3)  does  not  reach  the  brain 
before  the  transmitted  pressure  wave  from  pathways  1  and  2  completely  load  the  brain. 


^t  =  0ms  t  =  0.0136  ms  t  =  0.0271  ms  t  =  0.0541  ms 


t  =  0.06  ms 


Pressure  wave  from  Brain  region  not 
path  l  reaching  the  loaded  yet 

ventral  olfactory  bulb 


Pressure  wave  from 
path  3  reaching  the 
ventral  olfactory  bulb 


Wave  transmission 

Wave  transmission 
from  Path  I 


Interaction  of  pressure 


Figure  11:  Differential  loading  of  the  brain  caused  by  geometry  of  the  rat’s  head  and  animal  orientation  with 
respect  to  incident  shock  wave:  A.  Pressure  contour  plots  outside  the  rat  head,  B.  Corresponding  pressures 
inside  the  brain. 
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What  opportunities  for  training  and  professional  development  has  the  project  provided? 


Nothing  to  Report. 


How  were  the  results  disseminated  to  communities  of  interest? 


Nothing  to  Report. 


What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 


Tasks  in  the  Specific  Aim  2  (FY2)  include  determination  of  pathologies  associated  with 
primary  blast  exposure  corresponding  to  mild-to-moderate  TBI.  We  are  going  to  use  the 
survival  probability  dose  curve  developed  in  FY1,  to  identify  the  blast  parameters  associated 
with  50-75-100%  survival  for  detailed  pathological  and  proteomic  analysis. 

We  have  initiated  extensive  biochemical  characterization  of  various  brain  regions  extracted 
from  rats  exposed  to  blast  overpressure  in  the  130-180  kPa  range  and  sacrificed  24  hours  post 
injury  (per  experimental  design  included  in  the  proposal).  Goals  indicated  in  the  Statement  of 
Work  for  FY2  will  be  met  by  utilizing  qualitative  Western  blot  analysis  of  these  samples. 

This  analysis  will  be  supplemented  with  immunohistochemical  data  of  coronal  brain  sections 
stained  with  various  immunofluorescent  probes  to  identify  major  pathological  changes. 
Probing  will  be  performed  using  approximately  50  sections  per  brain  with  signal 
quantification  performed  in  the  entire  imaged  area.  For  this  task  Aperio  Versa  200  slide 
scanner.  Additional  information  regarding  the  extent  of  bTBI  will  be  obtained  from 
proteomics  analysis:  the  first  series  of  samples  is  currently  under  analysis  at  facilities  in  the 
Cancer  Center  at  Rutgers  NJMS. 
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4.  IMPACT 


What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 


Nothing  to  Report. 


What  was  the  impact  on  other  disciplines? 


Nothing  to  Report. 


What  was  the  impact  on  technology  transfer? 


Nothing  to  Report. 


What  was  the  impact  on  society  beyond  science  and  technology? 


Nothing  to  Report. 
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5.  CHANGES/PROBLEMS 


The  Project  Director/Principal  Investigator  (PD/PI)  is  reminded  that  the  recipient  organization  is 
required  to  obtain  prior  written  approval  from  the  awarding  agency  Grants  Officer  whenever 
there  are  significant  changes  in  the  project  or  its  direction.  If  not  previously  reported  in  writing, 
provide  the  following  additional  information  or  state,  “Nothing  to  Report,”  if  applicable: 


Changes  in  approach  and  reasons  for  change 


Nothing  to  Report. 


Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Describe  problems  or  delays  encountered  during  the  reporting  period  and  actions  or  plans  to 
resolve  them. 

Nothing  to  Report. 


Changes  that  had  a  significant  impact  on  expenditures 

Describe  changes  during  the  reporting  period  that  may  have  had  a  significant  impact  on 
expenditures,  for  example,  delays  in  hiring  staff  or  favorable  developments  that  enable  meeting 
objectives  at  less  cost  than  anticipated. 


Nothing  to  Report. 


Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents 

Describe  significant  deviations,  unexpected  outcomes,  or  changes  in  approved  protocols  for  the 
use  or  care  of  human  subjects,  vertebrate  animals,  biohazards,  and/or  select  agents  during  the 
reporting  period.  If  required,  were  these  changes  approved  by  the  applicable  institution 
committee  (or  equivalent )  and  reported  to  the  agency?  Also  specify  the  applicable  Institutional 
Review  Board/Institutional  Animal  Care  and  Use  Committee  approval  dates. 

Significant  changes  in  use  or  care  of  human  subjects 


Nothing  to  Report. 


~r~r 


Significant  changes  in  use  or  care  of  vertebrate  animals. 


Nothing  to  Report. 


Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  Report. 
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6.  PRODUCTS 


List  any  products  resulting  from  the  project  during  the  reporting  period.  If  there  is  nothing  to 
report  under  a  particular  item,  state  “Nothing  to  Report.” 

•  Publications,  conference  papers,  and  presentations 

Report  only  the  major  publication(s)  resulting  from  the  work  under  this  award. 

Journal  publications: 


1.  M.  Kuriakose,  M.  Skotak,  A.  Misistia,  S.  Kahali,  A.  Sundaramurthy,  N.  Chandra, 
Tailoring  the  Blast  Exposure  Conditions  in  the  Shock  Tube  for  Generating  Pure, 
Primary  Shock  Waves:  The  End  Plate  Facilitates  Elimination  of  Secondary  Loading  of 
the  Specimen,  PloS  one  11(9)  (2016)  e0161597. 

2.  V.  Mishra,  M.  Skotak,  H.  Schuetz,  A.  Heller,  J.  Haorah,  N.  Chandra,  Primary  blast 
causes  mild,  moderate,  severe  and  lethal  TBI  with  increasing  blast  overpressures: 
Experimental  rat  injury  model.  Scientific  Reports  6  (2016)  26992. 


Books  or  other  non-periodical,  one-time  publications. 


Nothing  to  report. 


Other  publications,  conference  papers,  and  presentations. 


Conference  presentations: 

1.  D.  Younger,  M.A.M.  Peringady,  D.  Haider,  N.  Prasad,  N.  Chandra,  Pathophysiological 
Changes  Due  To  Blast  Induced  Neurotrauma  Is  Effected  By  Animal  Orientation,  National 
Neurotrauma  Society  Meeting,  Lexington,  KY,  June  26-29,  2016  (J.  Neurotrauma,  2016, 
pp.  A87-A87). 

2.  B.  Swenson,  C.  Miao,  N.  Chandra,  B.  Pfister,  The  Spatial  And  Temporal  Deformation 
Pattern  Of  The  Brain  From  Blunt  Trauma,  National  Neurotrauma  Society  Meeting, 
Lexington,  KY,  June  26-29,  2016  (J.  Neurotrauma,  2016,  pp.  A22-A22) 

3.  S.  Kahali,  M.  Kuriakose,  M.  Skotak,  E.  Alay,  A.  Misistia,  N.  Chandra,  Artifact-Free 
Loading  Conditions  In  Compressed  Gas  Driven  Shock  Tube  With  End  Plate:  Focus  In 
Primary  Blast  Injury  Animal  Models,  National  Neurotrauma  Society  Meeting,  Lexington, 
KY,  June  26-29,  2016  (J.  Neurotrauma,  2016,  pp.  A93-A93)  * 

4.  G.  Ordek,  M.  Skotak,  N.  Chandra,  M.  Sahin,  Electrophysiological  Assessment  Of 
Blast-Induced  Injury  In  The  Rat  Cerebellum,  National  Neurotrauma  Society  Meeting, 
Lexington,  KY,  June  26-29,  2016  (J.  Neurotrauma,  2016,  pp.  A40-A41) 
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Website(s)  or  other  Internet  site(s) 

List  the  URL  for  any  Internet  site(s)  that  disseminates  the  results  of  the  research 
activities.  A  short  description  of  each  site  should  be  provided.  It  is  not  necessary  to 
include  the  publications  already  specified  above  in  this  section. 

Full  set  of  experimental  data  used  for  manuscript  preparation  is  provided  as  a  part  of 
publication  in  Plos  ONE:  doi:10.1371/journal.pone. 0161597. s009 


Technologies  or  techniques 

Identify  technologies  or  techniques  that  resulted  from  the  research  activities.  In  addition 
to  a  description  of  the  technologies  or  techniques,  describe  how  they  will  be  shared. 


Nothing  to  Report. 


Inventions,  patent  applications,  and/or  licenses 

Identify  inventions,  patent  applications  with  date,  and/or  licenses  that  have  resulted  from 
the  research.  State  whether  an  application  is  provisioned  or  non-provisional  and  indicate 
the  application  number.  Submission  of  this  information  as  part  of  an  interim  research 
performance  progress  report  is  not  a  substitute  for  any  other  invention  reporting 
reejuired  under  the  terms  and  conditions  of  an  award. 


Nothing  to  Report. 


Other  Products 


Theses: 

1.  A.  Misistia,  „Rodent  Testing  Device  Surrogate  for  Shockwave  Blast  Testing”,  2015. 

2.  S.  Kahali,  “Effect  of  endplate  on  the  blast  wave  profile  in  a  compressed  gas  shock  tube’ 
2015 


What  individuals  have  worked  on  the  project? 

Provide  the  following  information  for:  (1)  PDs/PIs;  and  (2)  each  person  who  has  worked  at  least 
one  person  month  per  year  on  the  project  during  the  reporting  period,  regardless  of  the  source 
of  compensation  (a  person  month  equcds  approximately  160  hours  of  effort).  If  information  is 
unchanged  from  a  previous  submission,  provide  the  name  only  and  indicate  “no  change.  ” 

Name:  Namas  Chandra 

Project  Role:  PI/PD 

Researcher  Identifier: 

Nearest  person  month  worked:  0.5 

Contribution  to  Project:  Design  of  the  research,  data  interpretation 

Name:  Maciej  Skotak 

Project  Role:  Research  Scientist 

Researcher  Identifier:  0000-0003-2584-7294 

Nearest  person  month  worked:  6 

Contribution  to  Project:  Exposure  of  animals,  sensor  implantation,  supervision  of  students 

&  staff,  data  analysis  (pressure  quantification),  report  preparation 

Name:  RamaRao  Venkata  Kakulavarapu 

Project  Role:  Research  Scientist 

Researcher  Identifier: 

Nearest  person  month  worked:  2 

Contribution  to  Project:  Overseeing  western  blot  and  histological  analysis,  supervision  of 

students  and  staff 

Name:  Eren  Alay 

Project  Role:  Laboratory  technician  (instrumentation) 

Researcher  Identifier: 

Nearest  person  month  worked:  4 

Contribution  to  Project:  Assistance  with  exposure  of  animals,  shock  wave  experimentation 

and  sensor  implantation,  sensor  calibration  and  instrumentation  maintenance 

Name:  Matthew  Kuriakose 

Project  Role:  Graduate  Student 

Researcher  Identifier: 

Nearest  person  month  worked:  4 

Contribution  to  Project:  Assistance  with  shock  wave  experimentation,  exposure  of  animals, 

biochemistry  work,  data  analysis,  manuscript  preparation 

Name:  Anthony  Misistia 

Project  Role:  Graduate  Student 

Researcher  Identifier: 

Nearest  person  month  worked:  4 

Contribution  to  Project:  Assistance  with  exposure  of  animals  and  sensor  implantation,  RTD 

development,  CAD  development,  data  analysis 
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Name:  Stephanie  Iring 

Project  Role:  Laboratory  technician  (biochemistry) 

Researcher  Identifier: 

Nearest  person  month  worked:  2 

Contribution  to  Project:  Western  blot,  staining  of  brain  sections 

Name:  Debanjan  Haidar 

Project  Role:  Graduate  student 

Researcher  Identifier: 

Nearest  person  month  worked:  2 

Contribution  to  Project:  Western  blot,  sectioning  of  brain  sections,  data  quantification 

Name:  Daniel  Younger 

Project  Role:  Graduate  student 

Researcher  Identifier: 

Nearest  person  month  worked:  1 

Contribution  to  Project:  Western  blot,  sectioning  of  brain  sections,  data  quantification 

Name:  Sudeepto  Kahali 

Project  Role:  Graduate  student 

Researcher  Identifier: 

Nearest  person  month  worked:  4 

Contribution  to  Project:  Numerical  simulations 


Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key  personnel 
since  the  last  reporting  period? 

If  there  is  nothing  significant  to  report  during  this  reporting  period,  state  “Nothing  to  Report.  ” 

If  the  active  support  has  changed  for  the  PD/PI(s)  or  senior/key  personnel,  then  describe  what 
the  change  has  been.  Changes  may  occur,  for  example,  if  a  previously  active  grant  has  closed 
and/or  if  a  previously  pending  grant  is  now  active.  Annotate  this  information  so  it  is  clear  what 
has  changed  from  the  previous  submission.  Submission  of  other  support  information  is  not 
necessary  for  pending  changes  or  for  changes  in  the  level  of  effort  for  active  support  reported 
previously.  The  awarding  agency  may  require  prior  written  approval  if  a  change  in  active  other 
support  significantly  impacts  the  effort  on  the  project  that  is  the  subject  of  the  project  report. 


Nothing  to  Report. 


What  other  organizations  were  involved  as  partners? 

Nothing  to  Report. 


24 


7.  SPECIAL  REPORTING  REQUIREMENTS 


COLLABORATIVE  AWARDS:  For  collaborative  awards,  independent  reports  are  required 
from  BOTH  the  Initiating  PI  and  the  Collaborating/Partnering  PI.  A  duplicative  report  is 
acceptable;  however,  tasks  shall  be  clearly  marked  with  the  responsible  PI  and  research  site.  A 
report  shall  be  submitted  to  https://ers.amedd.army.mil  for  each  unique  award. 

QUAD  CHARTS:  If  applicable,  the  Quad  Chart  (available  on  https://www.usamraa.army.mil) 
should  be  updated  and  submitted  with  attachments. 
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8.  APPENDICES 


a.  Reprint  of  the  paper  published  in  Scientific  Reports 


1 

www.nature.com/scientificreports 

SCIENTIFIC  REPpRTS 


open  Primary  blast  causes  mild, 

moderate,  severe  and  lethal TBI 
with  increasing  blast  overpressures: 
Experimental  rat  injury  model 

Published: 07  June  2016  V|kas  Wshrav  Ifeaej  Skotak1-',  Heather  Schuet2!,  Abi  HeHeH,  James  Haorah1  & 

Nam  as  Chandra1 


Injury  severity  in  blast  induced  Traumatic  Brain  Injury  <bTBl)  increases  with  blast  overpressure  (BOP) 
and  impulse  in  dose-dependent  manner.  Pure  primary  Mast  waves  were  simulated  in  compressed  gas 
shock-tubes  in  discrete  increments.  Present  work  demonstrates  24  hour  survival  of  rats  in  0-450  kPa 
(0-000  Pa-s  impulse)  range  at  10  discrete  levels  (60, 100, 130, 160, 190, 230, 250, 290, 350  and  420  kPa) 
and  determines  the  mortality  rate  as  a  non-linear  function  of  BOP.  Using  logistic  regression  model, 
predicted  mortality  rate  (PMR)  function  was  calculated,  and  used  to  establish  TBI  severities.  We 
determined  a  BOP  of  145  kPa  as  upper  mid  TBI  threshold  (5H  PMR).  Also  we  determined  140-220  kPa 
and  221-290  kPa  levels  as  moderate  and  severe  TBI  based  on  3544,  and  7044  PMR,  respectively,  while 
BOP  above  290  kPa  is  lethal.  Since  there  are  no  standards  for  animal  bTBl  injury  severity,  these 
thresholds  need  further  refinements  using  histopathology,  immunohistochemistry  and  behavior. 
Further,  we  specifically  investigated  mild  TBI  range  (0-145  kPa)  using  physiological  (heart  rate), 
pathological  (lung  injury),  immuno-hstochemical  (oxidative/nitrosative  and  Mood-brain  barrier 
markers)  as  well  as  blood  borne  biomarkers.  With  these  additional  data,  we  conclude  that  mid  bTBl 
occurs  in  rats  when  the  BOP  is  in  the  range  of  85-145  kPa. 


Exposure  to  blasts  is  one  of  the  leading  causes  of  trauma  experienced  by  military  personnel  as  a  result  of  wide¬ 
spread  use  of  high  explosives.  Our  unique  animal  models  of  primary  blast  waves  generated  in  compressed  gas 
shock  tubes  in  discrete  increments  shows  that  24-hour  survival  of  animals  depends  on  the  magnitude  of  blast 
overpressure  (BOP)u.  Blast-induced  traumatic  brain  injuries  (bTBl)  are  classified  as  primary,  secondary,  tertiary, 
and  quaternary3  *,  based  on  the  type  ofbiomechanical  loading  An  incident  pressure  of  a  shockwave  on  the  body 
within  the  time  duration  of  few  -to-  ten  milliseconds  causes  primary  blast  injury,  while  secondary  blast  injuries  are 
caused  by  impact  of  high  -  velocity  fragmentation  and  debris.  Tertiary  blast  injuria  result  when  the  body  is  vio¬ 
lently  accelerated  and  is  forced  to  impact  other  abiccts.  Quaternary  blast  injuries  are  caused  by  exposure  to  heat 
and  toxic  gases  released  resulting  from  explosive  detonation6.  The  present  study  is  focused  on  primary  bTBl  that 
was  recognized  as  a  separate  ncuropathological  condition  during  the  World  War  t  and  dubbed  as  ‘shell  shock'7. 
Soldiers  afflicted  with  ihdl  shock  in  the  battle  Ud  sustained  plethora  of  nriimlngjql  deficits  (or  even  death) 
without  any  visible  injuria  Long  after  the  shelling  had  ended1. 

The  emergence  of  bTBl  among  active  military  personnel  in  recent  Iraq  and  Afghanistan  war  gained  considera¬ 
ble  attention  and  remains  an  important  public  health  problem'*  B.  In  the  p>ast  decade,  the  Deportment  of  Defense 
reported  more  than  200.000  head  injuries  due  to  combat -related  incidents  and  in  non -deployed  environment11. 
The  severity  of  brain  injury  is  clinically  classified  as  mild11  “.moderate1*  severe*1*  a.  and  vegetative  state  TBI5 
as  per  15-pmdnt  Glascow  Coma  Scale  (GCS)  in  humans  (including  blast  TBI  cases14  M).  Over  150.000  of  these 
bead  injured  military  personnel  were  diagnosed  with  mild  Traumatic  Brain  Injury  (mTBI)  and  Post- Traumatic 
Stress  Disorder  (PTSD)  exhibiting  a  wide  range  of  neurological  and  psychological  symptoms*  V  Blast  mTBI  is 

!Center  for  Injury  Biomechanics,  Materials  and  Meticine  (OBM3),  Department  of  Biomedical  Engineering,  New 
Jersey  Institute  of  Technology,  Newark,  NJ  07102-1982,  USA. 7  Department  of  Pharmacology  and  Expen  mental 
Neuroscience,  university  of  Nebraska  Medea!  Center,  Omaha,  68196,  t£,USA.  These  authors  contributed  equally 
to  this  work.  Correspondence  and  recuests  for  materials  should  be  addressed  to  J  Jd.  (email:  james.haorah<£njrt-edu) 
orN-C  (email:  n  arras,  chandra®  nyt.edu) 


SCI ENT1FIC  REPORTS  1 6:26992  |  OCX:  10.103B/srsp26992 
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b.  Reprint  of  the  paper  accepted  for  publication  in  Plos  ONE 


(pjiPLOS  |  one 


RESEARCH  ARTICLE 

Tailoring  the  Blast  Exposure  Conditions  in  the 
Shock  Tube  for  Generating  Pure,  Primary 
Shock  Waves:  The  End  Plate  Facilitates 
Elimination  of  Secondary  Loading  of  the 
Specimen 

Matthew  Kurtakoee*,  Maclej  Skotak\  Anthony  MkUsha.  Sudeepto  K  shall, 

Aravlnd  Sundar  amurthy.  Names  Chan  A  a* 

Garter  tor topiy  Bormchanics.  Miuiik  and  Madana  (CCM3).  Oapa*nan(  d  Bomadcal  Engjnaanng. 

Nn w  Jesay  I raaiXe  d  Tachnotogy.  Namvk.  NJ.  071 02- 1982.  Unrtad  Stalaa  d  Armies 

w  Thoe*  autiarsoortatUad  aqualy  to  tiiswarfc. 

*  namaarhanckaOntLadu 


OPEN  ACCESS 


Abstract 


OtatarKultaai  U.  9o*k  M  Mtate  A 

S  SintaniTiiitty  A  Oianlra  S  ( 2D *)  Tafcrt-g  fra 
Efaat  Ej^xaur*  ConHkxv  in  fra  Shack  Tute  for 
Gaia teg  F\ji».  Prnury  Shark  Wms  Thi  End 
Ptatn  Facittaw  Elnfotaai  of  Scmtery  Lcad-g  of 
h«SpK*nai  PU>SCfE  np^aOtttSg? 

(tot©  O71^)urdponi01*iaB7 

Ed  tor  Firs  M  Kctensy  lEtenltyof  Rate 

iNirmswES 
teotmd  tey  17.2018 
Aceaatad:  Aigiat8. 2016 

Pubiahad  Sqtenfca  7. 20  tt 

Copyrl^it  0201;  KuiiAua  a  i  Thbkarcpar 
am  artkda  dtaWbufcd  urala  fri  bams  d  tia 
Oaalva Con irrrac  AtitxJcr  Uaa-aa.wHduwrria 

irratakte  is«  d-atufcr.  aid  lapradiclai  In  ary 
(radon  pnwtedfraatyraJaUtorinlsajRaara 

aadtal 

Qa  Ca  A  w  i  a  t*  Ity  S  titaaatf :  Dak  aa  potead 
mih  haSiffrteg  rfarraicr 

Furdrvj  Thk  wok  wassi^pcda:  by  grant  no 
VI® 900 1  anted  T>itnay  Bat  q  cry  Cdtate  tx 
AihalWunai  TBI  Modds  urirg  Fafc  Valcaac 
Shak  Tube*  rate  wd  (ten  fro  US .  Amy  tedtel 
Fteaach  ted  MdaW  Ctxranaxl 


The  end  plate  mounted  at  the  moutiofthe  shocktube  is  a  vereatite  and  effective  implement 
t>  contol  and  mrtigale  tie  end  effects.  We  have  performed  a  senes  erf  measurement  of 
incident  shoe*  wave  veto  cities  and  overpressures  foBowed  by  quanfficadion  of  impulse  val¬ 
ues  (integral  of  pressure  in  time  doman)  for  four  different  end  plate  configurations  (0.625, 2. 
4  inches,  and  an  open  end).  Shock  wave  characteristics  were  monitored  by  high  response 
rate  pressure  sensors  aBocated  in  six  positions  along  the  length  of  6  meters  long  229  mm 
square  cross  section  shodc  tube.  Tests  were  performed  at  tiree  shock  wave  in  tensities, 
which  was  controBed  by  varying  the  Mylar  membrane  fruckness  (0.02,0.04  and  0.06  inch). 
The  end  reflect)  r  plate  in  staBe  d  at  the  exit  of  tie  shock  tij  be  a  Bows  precise  control  over  tie 
intensity  of  retected  waves  penetrating  into  tie  shock  tibe.  At  tie  optimized  distance  of  tie 
lube  to  end  plate  gap  the  secondary  waves  were  entirely  elm  mated  from  tie  test  section, 
which  was  confirmed  by  pressure  sensor  at  T4  location.  This  is  pronounced  finding  for 
implementation  of  pure  primary  blast  wave  animal  model.  These  data  also  suggest  only 
deep  in  the  shoe*  tube  experimental  conditions  alow  exposure  to  a  single  shock  wave  tee 
of  artifacts  Our  restits  provide  detaied  insight  into  spa  tiotemporal  dynamics  of  shoe* 
waves  with  Friedander  waveform  generated  using  he  Bum  as  a  driver  gas  and  propagating 
in  the  at  inside  me  (turn  sized  tube .  Diffusion  of  driver  gas  (helium)  inside  tie  shock  tij  be 
was  responstile  for  velocity  increase  of  retected  shoe*  waves.  Numerical  simulations  com¬ 
bined  witi  experimental  data  suggest  the  shoe*  wave  attenuation  mechanism  is  simply  the 
expansion  of  tie  internal  pressure .  In  the  absence  of  any  o tier  pcetiiated  shockwave 
decay  mechanisms,  which  were  not  implemented  in  the  mod  el  the  agreement  between  the¬ 
ory  and  experimental  data  is  exceBent 


Cow  p«  tag  Iterate:  Thi  autxrshava  defer  to 

Eat  rc  oerpteg  hinds  «K 


PLOS  ONE  DOT:  10. 137  l^oumd-pona  .01  61 597  Saptamber7.2016 
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c.  Quad  chart 


Primary  Blast  Injury  Criteria  for  Animal/Human  TBI  Models  using  Field  Validated  Shock  Tubes 
14059001 

W81XWH-15-1-0303 

PI:  Namas  Chandra  Org:  New  Jersey  Institute  of  Technology  Award  Amount:  3,330,279 


Study/Product  Aim(s) 

•  Develop  two  master  dose-response  curves  correlating  mortality 
rate  with  blast  overpressure  (BOP,  D-450  kPa)  and  corresponding 
impulse  in  10  wk  old  Sprague  Dawley  rats 

•  Identify  and  characterize  the  mild  blast  TBI  thresholds  using 
developed  model  and  biochemistry:  characterize  oxidative  stress, 
BBB  damage  and  leakiness 

•  Evaluate  pathophysiology  of  rats  in  wide  range  of  BOPs  (0-4S0 
kPa) 

■  Develop  numerical  model  of  rat  brain  exposed  to  a  single  blast 
wave 


Approach 

Rats  were  exposed  to  a  single  blast  wave  at  NJIT  in  the  modular, 
multi-size  shock  tube  capable  of  reproducing  complex  shock  wave 
signatures.  We  evaluated  survival,  cardiac  and  pulmonary 
pathology  following  exposure  at  60, 100, 130, 160, 190, 230, 250, 
290,  350  and  420  kPa  peak  overpressure,  and  biochemical  assays 

in  60-130  kPa  range.  ..  .  _ 

Timeline  and  Cost 


Activities  CY 

16 

17 

18 

19 

Master  dose  response  curves 

Assess  mild-moderate  bTBI 

Effect  of  impulse  on  master  curves 

Establish  HIC  for  bTBI  (x-spedes) 

Estimated  Budget  ($K) 

$1,064 

$728 

$773 

$764 

Updated:  September  14th,  2016 


Goals/Milestones  (3  years  only) 

CY16  Goal:  Develop  Master  Dose  Response  Curves  for  10  w/o  SD  rats 
0  Evaluate  Mortality  &  Biomechanical  Loading  in  Wide  Range  of  BOP 
0  Determine  Biomechanical  Loading  of  Rat  Brain  in  Simulated  Blast 

□  Numerical  Simiiation  of  Brain  Injury 

CY17  Goal:  Assess  Pathologies  of  bTBI  24  hours  After  Exposure 

□  Evaluate  the  Lung  Injuries  Caused  by  Blast  Exposure 

□  Assess  Oxidative/Nitrosatrve  Stress.  BBB  Damage  and  Inflammation 
□Evaluate  Alterations  in  Brain  Proteome  After  Primary  Blast  Exposure 
CY18  Goal:  Effect  of  Blast  Impulse  on  Master  Dose-Response  Curve 

□  Establish  Master  Impulse  Dose-Response  Curve  at  Three  BOPs 

□  Protein  Expression  Due  to  Changes  in  Blast  Impulse 

□  Effect  of  Changes  in  Impulse  on  Loading  in  the  Rat  Brain 
C  omments/C  ha  llenges/  Iss  ues/Co  ncerns 

•  Timeline  change:  N/A. 

•  Spending  change:  N/A. 

Budget  Expenditure  to  Date 
Projected  Expenditure:  100% 

Actual  Expendrture:  75% 
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